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ABSTRACT 

Context. White dwarf masses in cataclysmic variables are difficult to determine accurately, but are fundamental for 
understanding binary system parameters, as well as binary evolution. 

Aims. We investigate the X-ray spectral properties of a sample of Intermediate Polars (IP) detected above 15keV to 
derive the masses of their accreting white dwarfs. 

Methods. We use data from the Swift/HAT instrument which during the first 2.5 yrs of operation has detected 22 
known intermediate polars. The X-ray spectra of these sources are used to estimate the mass of the white dwarfs. 
Results. We are able to produce a mass estimate for 22 out of 29 of the confirmed intermediate polars. Comparison 
with previous mass measurements shows good agreement. For GK Per, we were able to detect spectral changes due to 
the changes in the accretion rate. 

Conclusions. The Swift/~BAT detector with its combination of sensitivity and all-sky coverage provides an ideal tool to 
determine accurate white dwarf masses in intermediate polars. This method should be applied to other magnetic white 
dwarf binaries. 

Key words. X-Rays: binaries - stars: binaries: spectroscopic - stars: novae, cataclysmic variables 



1. Introduction 

Cataclysmic variables (CV) are binaries with a white dwarf 
(WD) accreting matter from a late main-sequence star 
which fills its Roche lobe (?). Depending on the strength 
of the magnetic field of the white dwarf, the matter flowing 
through the Lagrange point either forms a full or partial 
accretion disk or follows the magnetic field lines towards 
the magnetic poles of the white dwarf. Intermediate polars 
(IPs) are a sub-type of cataclysmic variables (see e.g. the 
review by ?) where a magnetic field of order B s=s 1-10 MG 
(?) disrupts the inner part of the accretion disk. From this 
disruption radius the accreting material follows the mag- 
netic field lines towards the magnetic poles and finally falls 
freely onto the surface of the WD, forming a strong shock 
near the surface of the WD (?). This so called post shock 
region is heated up to a high temperature, about 10-20 keV 
(?). The main cooling process is thermal Bremsstrahlung. 
The temperature of the post shock region and hence the 
spectral shape of the emitted X-ray radiation depend to 
first order only on the mass of the WD. This makes IPs 
especially interesting, because their X-ray spectrum allows 
a relatively simple estimate of the WD mass (see e.g. ???). 

The white dwarf mass is a fundamental parameter in 
cataclysmic binaries. It not only governs the dynamics of 
the orbital motion and the accretion flow, but also is a 
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main parameter characterizing the emission of the accretion 
region (see above). Moreover, it plays a fundamental role 
in the binary evolution, and accurate mass distributions of 
certain sub-classes can help in constraining the accretion 
history. 

Previous attempts to determine WD masses of IPs us- 
ing X-ray measurements have been controversial, since data 
from different X-ray instruments led to different mass esti- 
mates (? with Ginga data, or ? using RXTE/PCA, ASCA 
and Ginga data), and in all cases X-ray determined masses 
where substantially higher than masses derived from op- 
tical timing or spectroscopy. The classical example is XY 
Ari, an eclipsing IP (which is considered to be the most 
suitable object because the uncertainty in the inclination is 
very small): the RXTE spectrum yields 1.22 M Q , the ASCA 
data 1.27-1.40 M Q depending on GIS, SIS or a combination 
of the two detectors, and the Ginga data result in 1.3 M Q , 
while optical methods yield 0.78-1.03 M Q (?). 

Observations above 20 keV are crucial as they provide a 
clean signal not contaminated by other system components 
such as the accretion disc or the main-sequence star and 
it avoids complications due to cold or warm absorbers as 
well as fluorescence (line) emission. Moreover, at or above 
this energy, effects of absorption in the post-shock region 
are negligible. The only effect to take into account is reflec- 
tion which is known to be present in magnetic cataclysmic 
variables due to the detection of the 6.4 keV emission line. 
The only disadvantage was the weakness of the X-ray flux. 
Thus, a sensitive instrument must be used for this kind 
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Table 1. IPs detected in the BAT survey 3cr in the 15-55 keV band ordered by constellation. The first three columns 
give the BAT position and the offset to the optical position. The count rate was averaged over the entire 2.5 yr observing 
period. In the lower part of the table we report those sources which are not currently detected in the BAT survey; for 
those we give the optical position for which the upper limit was computed. 
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FO Aquarii 


334.481 
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12.4 ± 0.8 


16.6 


2.76 


1.99 


XY Arietis 


44.038 


19.441 


4.3 


5.3 ± 0.6 


7.9 
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1.55 


V405 Aurigae 
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6.7 


3.46 


1.60 
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1.77 
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PQ Geminomm 
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RX J2133.7+5107 


323.432 


51.124 


1.9 


10.9 ± 0.6 


18.2 


3.99 


1.42 


Swift J0732. 5-1331 


113.156 


-13.518 


2.8 


6.3 ± 0.6 


9.7 


2.61 


1.70 


IGR J00234+6141 


5.740 


61.686 


5.3 


2.0 ± 0.6 


3.0 


3.68 


1.72 


IGR J17303-0601 


262.590 


-5.993 


0.8 


14.0 ± 0.8 


12.7 


1.59 


2.78 


AE Aquarii 


310.04 


-0.87 




< 0.3 


0.4 


2.44 


2.05 


HT Camelopardaris 


119.26 


63.10 




< 0.1 


0.3 


3.52 


1.34 


DW Cancri 


119.72 


16.28 




< 1.5 


1.5 


2.09 


2.11 


UU Columbae 


78.05 


-32.69 




< 0.7 


1.7 


3.55 


1.50 


DQ Herculis 


271.88 


-45.86 




< 1.0 


2.5 


3.75 


1.56 


WX Pyxidis 


128.27 


-22.81 




< 1.4 


2.8 


3.16 


1.56 



1. Minimum flux to detect a source at 3cr at a given position in the 15-55 keV band. 



of observation. The Swift Burst Alert Telescope (BAT, ?), 
with an unprecedented sensitivity in the 15-200 keV range, 
opened a new window for the study of the hard X-ray sky, 
combining high sensitivity and all-sky coverage. 

The present paper is organized as follows: in section 2 
we introduce the observations. The results of the measure- 
ment can be found in section 3 which is divided in three 
subsections to discuss first the spectral analysis, second the 
accretion-rate dependence observed in GK Per, and finally 
the comparison with previous methods of mass determina- 
tion (both with X-ray and non X-ray methods). The con- 
clusions are given in section 4. 

2. Data and models 

2.1. The Sw\h/BAT X-ray survey 

BAT is a coded mask telescope with a wide field of view 
aperture (FOV, 120° x 90° partially coded) sensitive in the 
energy range of 15-195 keV. BAT continuously monitors the 
hard X-ray sky surveying gamma-ray bursts. Results of the 
BAT survey (??) show that BAT reaches a sensitivity of 
«1 mCrab in a 1 Ms exposure. With its sensitivity and 
the large amount of exposure already accumulated over the 
whole sky, BAT is a very good instrument for studying IPs 
above 20 keV. 



In the analysis presented here we use BAT data to derive 
spectra of known IPs. This was done by cross-correlating 
our list of source candidates with catalogs of IPs. 

For the source detection step, we used all the available 
BAT data taken from January 2005 to March 2007. An en- 
ergy range of 15-55 keV was chosen. The lower limit is given 
by the energy threshold of the detectors. The upper limit 
is taken to avoid the presence of strong background lines 
which could worsen the overall sensitivity. The data screen- 
ing was performed analogously to the method presented in 
?. The final all-sky image is obtained as the weighted av- 
erage of all the shorter observations. The average exposure 
time in our image is 3 Ms, with 1.3 Ms minimum and 5 Ms 
maximum exposure times. The final image shows Gaussian 
normal noise and we identified source candidates as excesses 
above the 3er level. All the candidates are then fit with 
the BAT point spread function (using the standard BAT 
tool batcelldetect) to derive the best source position. The 
BAT position error is between 3-7 arcmin, depending on 
the source intensity (see ?). The identification of IPs from 
this detection list was done using the online catalog main- 
tained by K. Muka{3 and allowing for each source's position 
error. 



1 http:/ /asd. gsfc.nasa.gov/Koji.Mukai/iphome/iphome. html, 
version of March 2007 
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Fig. 1. Swift/BAT spectrum of GK Per, averaged over 2.5 
yrs, fitted with the model described in ?. 



Out of the 29 known IPs to date, Swift/BAT detects 22 
objects (see Table [J). The sources that fall below the BAT 
sensitivity (AE Aqr, HT Cam, DW Cnc, UU Col, DQ Her 
and WX Pyx) are likely either very far away, or are in a 
very low accretion state. All six sources except HT Cam 
(RXTE detection; ?) were also not detected in RXTE or 
Integral/IBIS observations. 

For each of the detected IPs, we extracted a 15-195 keV 
spectrum with the method described in ?. Here we recall 
the main steps. For a given source, we extract a spectrum 
from each observation for which the source is in the field of 
view. These spectra are corrected for residual background 
contamination and for vignetting. The per-pointing spectra 
are then (weighted) averaged to produce the final source 
spectrum. Thus, the final spectrum represents the average 
source emission over the entire survey period (2.5 years) 
although individual coverage intervals may differ. 

We also extracted the light curves of all detected IPs. 
In most of the cases, the sources are too faint to assess any 
variability. The only exception is GK Per (see Fig. [2]) and 
this source will be discussed in detail in Section P 



this model does not include reflection, ? have shown that 
including the reflection changes the fit parameters and the 
corresponding white dwarf masses only slightly. The main 
difference between the model of ? and that of Aizu (1973) is 
that ? uses the cooling functions of ?, which takes thermal 
Bremsstrahlung, lines and b-f transitions into account. 

In the following we have used the model of ?. The free 
parameters are the mass of the white dwarf and a normal- 
ization constant. Since we fit spectra only above 15 keV, 
we do not apply any foreground absorption. 

3. Results 

3.1. Spectral analysis 

In order to derive a WD mass estimate we fitted the model 
of ? to the BAT IPs using XSPEC. An example is given in 
Figure [TJ The detection of the IP emission up to ^100 keV 
allows an accurate estimate of the bremsstrahlung temper- 
ature. Since the X-ray spectrum depends primarily on the 
WD mass, using the model of ? we are able to determine 
directly an estimate (and its uncertainty) of the WD mass. 

With the exception of MU Cam, all spectral fits show 
acceptable reduced x 2 ( see Table [2] which include the er- 
rors at the 90% confidence level; the number of degrees 
of freedom is always 10, based on 12 channels and 2 fit 
parameters). We refrained from attempting to fit more so- 
phisticated models. Note that, in all the cases, a power law 
does not produce a good fit because of the clear curvature 
of the data. 



3.2. Accretion- rate dependence of GK Per 

In the simplistic case, it is assumed that the accreted mass 
falls freely from infinity onto the white dwarf. In reality, 
however, the matter does not come from infinity, but falls 
from the Alfven radius Pm, which in turn depends on the 
mass accretion rate (from ?): 



Ra 
Rwd 



2.3 



M 



-2/7 



10 20 g/s 



WD 



-1/7 



%D \ 

10 9 cm / 



5/7 



( B 



Vl0 6 G 



2.2. Spectral model 

The main cooling mechanism in the post-shock region is 
through bremsstrahlung, thus the simplest method is to 
fit a Bremsstrahlung spectrum to the data. The struc- 
ture of the post-shock region was investigated first by 
?, who provided a conversion formula to determine the 
white dwarf mass based on the effective temperature of 
the Bremsstrahlung spectrum. However, in order to esti- 
mate the mass of the accreting white dwarf more accu- 
rately, one needs a more careful modeling of the temper- 
ature and emissivity distribution in the post-shock region. 
? took into account cyclotron cooling, which can be im- 
portant for polars, but was shown to be unimportant for 
objects with surface magnetic fields less than 10 MG (like 
IPs). ? derived a simple analytical model under the as- 
sumption of constant pressure in the post-shock region. ? 
extended this model to take into account how the pressure 
grows towards the white dwarf surface. The emergent spec- 
trum thus is the sum of many local Bremsstrahlung spectra 
with their appropriate temperatures and pressures. While 



Thus, including this dependency, the white dwarf mass 



is: 



Mwd(^a) = M WD (oo) 1 



Rwd \ 

~rTJ 



(2) 



Therefore, the higher the accretion rate, the smaller is 
Ra, and thus the smaller the mass of the white dwarf which 
is deduced under the assumption of infinite origin. 

There is yet another dependency which makes a state- 
ment on the deduced mass relative to the real mass some- 
what uncertain. Accretion is only possible if Ra < Rn (?), 
where Rn is the "corotation" radius, i.e. the radius where 
the Keplerian angular velocity is equal to the rotation ve- 
locity of the WD (otherwise the mass would be thrown out 
of the system): 



Rn 
Rwd 



= 2.3P p 2/3 (min) 



A/- 



WD 



MfT 



1/3 



Rwd 
10 9 cm 



(3) 



Thus, fast rotating systems require higher accretion 
rates at the same magnetic field strength, and we measure 
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Fig. 2. Light curves of GK Per in the 15-195 keV band with a binning of 5 days. The middle panel shows the full 
observation while the left panel shows the outburst at the beginning and the right panel that at the end of the covered 
observing period. 



Table 2. Detected IPs with masses from XSPEC fit with the model from ?. The temperature is from a bremsstrahlung 
fit, and is provided for comparison. 



Name kT Mass Fi4_ig 5 Distance' 4 Li4_igs M x?cd 

(keV) (M ) (10~ n erg/cm 2 /s) (pc) (10 32 erg/s) (10 16 g/s) 



FO Aqr 


15.2±1.6 


0.61 


± 


0.05 


4.7 


± 


0.5 


400 1 


9.00 


29.0 


1.4 


XY Ari 


29.1±6.2 


0.96 


± 


0.12 


2.3 


± 


0.4 


270 2 


2.01 


2.8 


1.3 


V405 Aur 


22.1±4.2 


0.89 


± 


0.13° 


2.2 


± 


0.5 








0.5 


MU Cam 


19.6±4.7 


0.74 


± 


0.13 


1.3 


± 


0.3 








1.7 


V709 Cas 


30.0±2.5 


0.96 


± 


0.05° 


7.2 


± 


0.6 








1.0 


V1025 Cen 


10.8±7.7 


0.46 


± 


0.31 


0.4 


± 


0.3 








0.6 


BG CMi 


18.9±5.5 


0.67 


± 


0.19 


2.0 


± 


0.4 


700 3 


11.73 


32.3 


0.7 


TV Col 


21.6±2.4 


0.78 


± 


0.06 


5.0 


± 


0.5 


386 4 


8.92 


18.8 


0.9 


TX Col 


17.0±3.2 


0.67 


± 


0.10 


1.7 


± 


0.3 


550 5 


6.15 


16.9 


0.7 


V2306 Cyg 


20.5±5.7 


0.77 


± 


0.16 


1.2 


± 


0.3 








1.0 


DO Dra 


12.0±2.6 


0.50 


± 


0.11° 


1.0 


± 


0.3 


155 6 


0.29 


1.3 


0.7 


PQ Gem 


16.8±2.8 


0.65 


± 


0.09 


2.5 


± 


0.4 








1.0 


EX Hya 


19.4±4.4 


0.66 


± 


0.17 


2.1 


± 


0.5 


64.5 7 


0.10 


0.28 


0.4 


NY Lup 


39.1±4.3 


1.09 


± 


0.07 


7.1 


± 


0.6 


690 s 


40.45 


25.3 


1.1 


V2400 Oph 


21.1±2.8 


0.81 


± 


0.10 


4.3 


± 


0.6 








1.0 


GK Per 


26.2±5.4 


0.90 


± 


0.12 s,c 


2.8 


± 


0.5 


340 9 


10.24 


16.2 


0.8 


AO Psc 


13.1±1.9 


0.55 


± 


0.06 


2.8 


± 


0.4 


420 10 


10.13 


38.6 


0.7 


V1223 Sgr 


17.0±1.3 


0.65 


± 


0.04 


8.5 


± 


0.8 


527 11 


28.25 


82.0 


1.0 


V1062 Tau 


18.3±5.8 


0.72 


± 


0.17 


1.6 


± 


0.4 


500 12 


4.79 


11.6 


0.4 


RX J2133 


27.1±2.9 


0.91 


± 


0.06 


4.6 


± 


0.4 


>600 13 






0.6 


Swift J0713 


21.7±4.1 


0.79 


± 


0.11 


2.5 


± 


0.4 








0.9 


IGR J00234 


22.8i??° 


0.85 


± 


0.39 


0.7 


± 


0.3 








0.8 


ICR J17303 


37.1±4.4 


1.08 


± 


0.07 


6.6 


± 


0.6 








1.1 



A References to the distances: 1. ?; 2. ?; 3. ?; 4. ?; 5. ?; 6. ?; 7. ?; 8. ?; 9. ?; 10. ?; 11. ?; 12. ?; 13 ?. 
B All numbers correspond to the measured spectrum during the low-flux period; but see section 3.2. 
c The derived mass is likely a lower limit due to the fast rotation of the WD; see section 3.1. 



smaller masses than they have in reality. Because of this, 
we think we measure only a lower limit for the mass of the 
following four systems (with their WD rotation period in 
parenthesis): GK Per (5.85 min), V709 Cas (5.12 min), DO 
Dra (8.83 min), V405 Aur (9.09 min). 

Fig. [2] reveals that GK Per had an outburst at the be- 
ginning (January 2005) as well as at the end (Dec. 2006 
- >Mar 2007) of the covered observing period. GK Per 
is well known for its periodic outbursts, which take place 
every three or four years (??). ? mention an increased opti- 
cal brightness in Dec. 2006, about 0.6 mag above its usual 
quiescent level, which is consistent with our detection of in- 
creased X-ray emission. Since GK Per is among the bright- 



est IPs in our sample, we are able to split the data into 
three different periods, namely the low- flux period as well 
as the two outbursts, and derive the mass separately for 
these periods. 

The results in Table [3] show that we measure lower 
masses at higher flux. While our data are not sufficient to 
demonstrate that the mass estimates quantitatively follow 
that expected on the basis of the expected change in R^, 
the change is suggestive of being due the dependence on 
the accretion rate according to equations [1] and [U Note the 
systematically higher mass estimate if we were not been 
able to distinguish the different flux levels. Only the mass 
derived during the low-flux (low-accretion) period of GK 
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Table 3. Comparison of mass estimates for different sec- 
tions of the light curve of GK Per. 



Part of light curve 


Mass (M Q ) 


Xrcd 


Low-flux period, without outbursts 


0.90±0.12 


0.78 


Jan. 2005 outburst 


0.74±0.05 


1.66 


Mar. 2007 outburst 


0.67±0.03 


1.24 


Sum of the total period 


0.78±0.05 


0.79 



M =1.1 

IVI WD 



Per should be considered - the last row in Table [3] is for 
illustrative purposes only. The values as derived during the 
outbursts are similar to those of ?, suggesting that those 
data were also taken during periods of increased accretion. 
We finally note that for all other sources of TableOthe mass 
estimates are unlikely to be affected by high accretion rates, 
i.e. the source being in outburst. 

Assuming that the magnetic field (and the mass) does 
not change significantly, we attempted to determine the 
accretion rate via equations [TJ [5] and using the rotation 
period of 351 sec ?. In Fig. [3] the three mass estimates 
from Tableware shown together with the generally adopted 
values for the mass accretion rate of GK Per in its different 
states (?), namely 4xl0 16 g/s at quiescence and 23xl0 16 
g/s in outburst. We also add the measurement of M = 0.6 
Mq from ?. In order to be consistent, all four measurements 
should fall onto one particular line for the magnetic field of 
the white dwarf. Testing white dwarf masses between 0.6- 
1.3 Mq, we find that the mass of the white dwarf in GK 
Per must be definitely higher than 0.9 M©. A mass of 1.15 
Mq provides the best match, and we present the graphs for 
1.1 and 1.2 M© to show how sensitive the mass estimate 
is (note, however, that this hinges on the assumption of 
the accretion rate). For 1.1 Mq the data points are flatter 
than the corresponding curve of B = 0.2 MG, while for 1.2 
Mq the data points are already steeper than the nearest 
B=const. curve. The online material provides the remaining 
figures with white dwarf masses of 0.9 Mq, 1.0 M Q and 1.3 
Mq . Our accretion rate corrected mass of the white dwarf 
of GK Per of 1.15 M Q compares well with the lower limit of 
0.9 M derived by Morales-Rueda et al. (2002). However, 
the implied magnetic field of 0.2 MG is at the low end of 
what is believed to be the canonical magnetic field strength 
of IPs, namely 1-10 MG (?). 

The canonical way of deriving the accretion rate is 
via the X-ray luminosity, i.e. using the relation Lx = 
0.1 x e x G x Mwd x M/Rwd, where e is the bolometric 
correction to the X-ray luminosity in the 14-195 keV band, 
the factor 0.1 accounts for the ratio of X-ray to total ac- 
cretion luminosity, and we assume the conversion efficiency 
of gravitational binding energy to radiation being equal to 
1 for WD surface accretion, as well as isotropic emission. 
For a bremsstrahlung spectrum of 20 (40) keV, e = 0.33 
(0.53). Using our mass estimate of 0.90 (1.15) M Q and the 
corresponding radius of 6200 (3500) km (according to the 
relativistic Fermi gas), we obtain M— 1.7 (0.7) xlO 16 g/s. 
This is somewhat lower than according to the typical val- 
ues used above, and, if true, would suggest an even lower 
magnetic field. It seems obvious that the combination of 
the accretion rate estimate from Lx with the dependencies 
according to eqs. 1-3 is not fully consistent with our general 
belief of B ~ 1-10 MG in IPs. 
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Fig. 3. Relation between measured mass and accretion 
rate for different magnetic fields for Mvfd=1-1 M q and 
Mwd =1-2 M Q . The values of this work are shown (filled 
circles) with accretion rates from ?. The parallel curves are 
for different values of the magnetic field, namely 0.1, 0.2, 
0.3, 0.4, 0.6, 0.8, 1, 2, 3, 5, 7 and 10 MG, respectively (from 
bottom to top). The open square represents the measured 
mass and its corresponding accretion rate derived by ?. The 
accretion rate has its first limit from the constraint that the 
system must have an accretion rate lower than 0.01T0 20 g/s 
at quiescence, because of its life time. 



3.3. Comparison with other methods of mass determination 

The canonical method of determining the WD mass in CVs 
is by measuring the radial velocities of optical emission 
lines. However, this method suffers from two main prob- 
lems. First, the uncertainty in the inclination of the system 
leads to a wide range in mass. Second, in magnetic systems 
the velocities along the collimating field lines dominate the 
radial velocities. 

An alternative method is to measure the mass ratio and 
infer the mass of the secondary by modeling the infrared 
light curves (e.g. ?). Again, the uncertainty in the inclina- 
tion limits the accuracy of the mass determination. 
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M WD from X-ray (this work), M 



Fig. 4. Comparison of the measured white dwarf masses 
with the ones from other methods (values from Table Q| . 
The source names are abbreviated; the labels of V1223 Sgr 
and EX Hya overlap. 



Eclipsing systems offer a unique opportunity to measure 
the mass more accurately due to the relatively small uncer- 
tainty in the inclination of the binary system. Moreover, in 
eclipsing systems the radius of the white dwarf can be deter- 
mined, thus providing yet another method (though rarely 
applicable due to the few eclipsing systems) to determine 
the WD mass. 

3.3.1. The whole sample 

In Table 2] and Fig. [J] we compare the mass estimates ob- 
tained by our analysis to previous mass estimates via radial 
velocity studies. In general, there is good agreement, with 
no systematic effect visible. A detailed comparison with re- 
spect to individual sources is given below. 

While we compare our results to those obtained with 
IBIS onboard INTEGRAL (see below), we choose not to 
discuss the differences (Table [4} of our results to those ob- 
tained by Ginga, ASCA and RXTE due to the systematic, 
and unexplained, differences for XY Ari (?). In the next 
section we discuss each object. 

3.3.2. Discussion of individual objects 

The mass of 0.96±0.12 M© for XY Ari matches very well 
the value of 1.04±0.13 M Q from ?. This is particularly en- 
couraging given the systematic differences in previous esti- 
mates as discussed in ? - see the Introduction. Since XY 
Ari is the only eclipsing IP, its optically determined mass 
is generally assumed to be secure. The fact that our mass 
estimate matches so well suggests that our analysis does 
not suffer from major systematic uncertainties. 

? derives a mass range of 0.6-1.0 M Q for the WD in 
BG CMi. Our value of 0.72±0.14 M Q falls in that range. 

For the system TV Col, ? determined a mass of 
0.75±0.15 Mq, which agrees very well with our estimate 
of 0.78i0.06 M(D. 



In DO Dra our WD mass estimate is 0.50±0.11 M . ? 
find a mass of 0.83±0.10 M Q , which is substantially differ- 
ent to our value. ? also find a mass of 0.83±0.10 M©. The 
difference could be the result of the short spin period of the 
WD (see subsection l3.2p . thus making our estimate a lower 
limit. However, our mass estimates of the other WDs with 
short spin period are not as far below the optical estimates. 

Using data of the first months of BAT observations, a 
significantly higher mass (0.83±0.17 M Q ) is deduced for 
EX Hya than for the rest of the time - although the count 
rate does not differ in these periods. If the first 100 days 
are ignored, a mass of 0.66±0.13 M Q is derived. While we 
cannot completely exclude a statistical fluctuation which 
leads to a change in the spectrum, one possible physical 
reason for the strange behavior could be that this system 
has a much lower mass accretion rate than the others. In 
the model of ? cyclotron cooling plays a role for a < 0.01 
g/s/cm 2 , so if the accretion rate were lower, this could pro- 
vide an explanation. Hoogerwerf et al. (2004) measured a 
mass for the WD in EX Hya of 0.49±0.13 M , which is 
comparable within the errors with our value. 

In the system NY Lup we obtain a mass of 1.12±0.06 
Mq. This value is consistent with the lower limit of 0.5 M Q , 
which Martino et al. (2006) determined. 

For V1223 Sgr a WD mass of 0.65±0.04 M was de- 
rived. ? gives a range of 0.4-0.6 M Q , which fits our value. 

The accretion-corrected mass of 1.15±0.15 M Q for GK 
Per is within the error bars of the values from other meth- 
ods: 0.9±0.2 Mq from Crampton et al. (1986) and the lower 
limit M>0.87±0.24 M Q from ?. 

3.3.3. Comparison with INTEGRAL/IBIS data 

? and very recently ? provided a table with temper- 
atures derived from Integral/IBIS data. The method 
of Aizu (1973) provides a direct conversion from the 
bremsstrahlung temperature to mass of the white dwarf, 
which we use to derive the masses as listed in Table [5] 
The comparison of these masses to those derived here using 
Swift/BAT shows, in general, very good agreement (Table 
[5] The only outlier is MU Cam where our temperature (and 
mass) is substantially greater. However, in a recent analy- 
sis of XMM-Newton and INTEGRAL data, ? conclude that 
the estimate of ? seems too low. These authors attempted 
a combined XMM-INTEGRAL fit, and arrive at a temper- 
ature of 35±10 keV, which they describe as "convincing 
agreement with the INTEGRAL data points" . Though the 
Swift/BAT spectrum for MU Cam is the worst in terms of 
reduced x 2 > we are confident that our measurement does 
not suffer from large systematic uncertainties. 

4. Conclusions 

An accurate determination of white dwarf masses in IPs 
is important for our understanding of the accretion geom- 
etry and binary evolution. However, this has so far been 
complex because of the difficulty in measuring the radial 
velocities and inclination of the systems. The main advan- 
tage of using X-ray measurements is that they allow a direct 
measurement of the WD mass without relying on the above 
parameters. 

Indeed, since bremsstrahlung is the main cooling mech- 
anism in the post-shock region, it is possible to relate the 
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Name 


BAT 1 


RXTE 2 


RXTE 3 


GINGA 4 


ASCA 5 


Other methods 6 


References 7 




M 


M 


M 


M Q 


M 


M 




FO Aqr 


0.61±0.05 


0.60±0.05 


0.88 


0.92 


1.05 






XY Ari 


0.96±0.12 




0.97 






1.04±0.13 


[1] 


V405 Aur 


0.89±0.13 


0.90±0.10 


1.10 




>0.54 




V709 Cas 


0.96±0.05 


0.90±0.10 


1.08 










V1025 Cen 


0.46±0.31 




0.60 










BG CMi 


0.67±0.19 


0.85±0.12 


1.20 


1.09 




0.8±0.2 


[2] 


TV Col 


0.78±0.06 


0.84±0.06 


0.96 


1.30 


0.51 


0.75±0.15 


[3] 


TX Col 


0.67±0.10 


0.70±0.30 


0.73 


0.48 


0.66 




DO Dra 


0.50±0.11 


0.75±0.05 








0.83±0.10 


[4] 


PO Opm 

1 \c£ VjrcIIl 








1.29 






EX Hya 


0.66±0.13 


0.50±0.05 


0.45 


0.46 


0.48 


0.49±0.13 


[5] 


NY Lup 


1.09±0.07 










>0.5 


[6] 


V2400 Oph 


0.81±0.10 


0.59±0.05 


0.71 




0.68 






GK Per 


0.90±0.12 


0.59±0.05 






0.52 


0.9±0.2, >0.87 


[7], [8] 


AO Psc 


0.55±0.06 


0.65±0.05 


0.60 


0.56 


0.40 




V1223 Sgr 


0.65±0.04 


0.95±0.05 


1.10 




1.28 


0.4-0.6 


[2] 


V1062 Tau 


0.72±0.17 


1.00±0.20 


0.90 









References: 1: this work, BAT; 2: ?; 3: ?; 4: ?; 5: Ezucka,Ishida (1999); 6: Other methods of determination of mass and the last 
column are the references of them. 7: References: [1] ?, [2] ?, [3] ?, [4] ?, [5] ?, [6] ?, [7] ?, [8] ?. 

Table 5. Comparison of masses and effective temperatures with a remsstrahlung model from BAT (this work, model: ?) 
and INTEGRAL/IBIS by (?) and (?). 1: Masses, derived from the temperatures with the model of ?. 2: Bremsstrahlung 
temperature. 



Name BAT/this work IBIS/Barlow IBIS/Landi 

mass (M ) kT (keV) mass (Me) 1 kT (keV) 2 mass (Me) 1 kT (keV) 2 



FO Aqr 


0.61±0.05 


15.2±i;S 






1 06+ ' 38 

1 - uo -0.51 


29 7+ 70 - 1 

' -16.6 


MU Cam 


0.74±0.13 


19.6±4.7 


0.37±0.17 


8.1±4.7 






V709 Cas 


0.96±0.05 


30.0±2.5 


0.88±0.06 


23.3±2.2 


n qr+0.07 

u.ao_ 08 


25.6l 2 ;I 


NY Lup 


1.12±0.06 


39.1±4.3 


0.99±0.06 


27.1±2.2 






V2400 Oph 


0.81±0.10 


21.1±2.8 


0.73±0.05 


18.6±1.4 






GK Per 


0.90±0.12 


26.2±5.4 


1.00±0.40 


28.7±15.6 


0.87l°; 2 ? 


23.0t 9 6 . 2 


V1223 Sgr 


0.65±0.04 


17.0±1.3 


0.74±0.07 


18.8±1.2 






RX J2133 


0.91±0.06 


27.1±2.9 


0.90±0.12 


23.8±4.3 






ICR J00234+6141 


0.85±0.39 


22.81??'° 


0.64±0.18 


15.9±5.1 






IGR J17303-0601 


1.08±0.07 


37.1±4.4 


0.98±0.14 


26.7±4.8 




31.6tft 7 



temperature of the IP's X-ray spectrum to the WD mass 
(e.g. ??). Using this method, we determined WD masses 
for a complete sample of 22 IPs detected in the ongoing 
Swift/BAT survey. For 6 objects, this represents the only 
available mass estimate. 

All mass estimates are in good agreement with previ- 
ously estimated WD masses and in general, thanks to the 
good signal-to-noise ratio, are more accurate. In particular, 
we reach surprisingly good agreement for the only eclips- 
ing IP, XY Ari, suggesting that our analysis does not suffer 
from major systematic errors. 

Detailed analysis of the light curve of GK Per shows 
that BAT detects the dependence of the source flux on the 
accretion rate. While this has been expected for a long time, 
to our knowledge this effect has never yet been supported 
by observational evidence. 
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